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ABSTRACT 


In this work, thermal fatigue analysis of pressure 
vessels with hemispherical, elliptical, torispherical 
and toriconical heads has been studied using finite ele- 
ment method. Inside of the vessel is subjected to pre- 
ssure and cylic temperature variation of heating, hold- 
ing and cooling; and outside is insulated. Weight of 
the vessel has been considered and isoparametric eight 

noded quadrilateral ring elements have been used. 

* 

Stress variation is studied in all the vessels. >.as 
function of heating/cooling rate and hold time. It is 
observed that maximum stress variation occurs in the 
skirt region near the head and is tie dominant factor 
for the thermal fatigue life of the vessels with all 
heads. Vessel with hemispherical head is the best design, 
and others in this order are elliptical, torispherical and 
toriconical heads, 

Finally the design changes axe considered to reduce 
the stress variation in the skirt, Qne is to change the 
position of the skirt and other is to provide a portion 
of it a temperature variation as that of the inside of 
the vessel. The second .method reduces the stress variation 
considerably. 



OKAPIS'! I 
I IT ISO LUG TI OH 

In recent years of research and development considerable 
attention has been focused on the design of pressure vessels 
because these are the means by which energy systems controlled, 
chemical and petroleum processes operated. Modern chemical 
processing plants and power generation industry demand the use 
of very high pressures and temperatures. Many chemical pro- 
cessing pressure vessels are subjected to a large number of ; 

i 

temperature cycles, hence to stress cycles. Ultimately the j 

j 

vessel is liable to fail 'by thermal fatigue. Thus in such 

i 

cases the designer is charged with the task of doing thermal 
fatigue analysis. This can only be done by a careful stress j 
analysis of the entire vessel. 

1,1 PREVIOUS WORK 

Pressure vessels and its components are designed based 
on the standard codes like ASKE, US §tc. However, the design- 
er is often faced with complex loadings which are not covered 
by the cades. Thus certain designs fall beyond the scope of 
the code. One such loading is that of the inside surface of 
the vessel being subjected to cyclic temperature variation. 

This will induce transient thermal stresses which depend on 
the temperature attained and the thermal gradients caused. It 
is not possible to estimate the temperatures and the corres- 
ponding thermal, stresses, as closed form relations are not 
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available. Go a design of the vessel is completed on the basis 
of mechanical loads. 

When once this design is completed, a detailed thermal 
stress analysis is needed when the vessel is subjected to 
thermal transients. In thermal fatigue problems, the designer 
is interested in knowing what is the critical section, where the 
stresses are going to be maximum. Also for calculating the 
allowable number of thermal cycles, knowledge of maximum alter- 
nating stress or strain amplitude is essential. 

Stress analysis of pressure vessels can be performed by 
the following methods: 

1 . Analytical Methods 

2. Experimental Methods which include 

a. Strain gauge measurement 

b. Photo elastic method 

c. Hr it tie coating method 

3. numerical Techniques. 

Analytical methods require a general solution of the equi- 
librium equations. But sometimes the problem becomes too complex' 
and beyond analytical solution. Then one has to resolve to experi- 
mental methods. But the experimental methods pose a lot of diffi- 
culties in the fabrication of the model and creating the required 

of 

type of loading. This needs a lot/skill and is highly costly in 
time and money. 

Thus in recent times, numerical Techniques such as finite 
element method has been used for the stress analysis* In EEM, 
the continuum is divided into a finite number of elements* A 
solution for the displacements is assumed initially in the form 
of a polynomial without giving any consideration to the boundary 
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conditions of the problem* ITov; -using either variational princi- 
ples or weighted residual methods, typical finite element equa- 
tions between stiffness, displacements and load matrices are 
formed. (These elemental equations are suitably assembled and 
global equations are formed. Then by applying the boundary 
conditions, the system of equations .are solved for the displace- 
ments and loter strains and stresses are found at various points. 
In this way a general PPM program can handle a variety of pro- 
blems of Solid Mechanics, Heat Transfer and Pluid Mechanics* 

The development of modern computers with tremendous speeds of 
calculation made PEEK to emerge as a powerful numerical technique 
in the field of stress analysis* The results from the PEM are in 
excellent agreement with theoretical exact solutions available 
for some problems. 

larger £lj reported tb i 1 since most pressure vessels are 
subjected to limited number of pressure and temperature cycles 
during the life time, considerable design effort could be saved 
by defining the conditions which do or do not require a fatigue 
evaluation. The author in his studs'- concluded that a fatigue 
evaluation is needed if during normal operation, 

1. Pressure fluctuates through a range exceeding 
twenty-five percent of the design pressure, or 

2. The temperature difference between any two adjacent 
points in the vessel changes by more than 22°G (40®P). 

This mokes it possible to eliminate a great deal of detailed 
calculations for vessels, that are subjected to pressure and tem- 
perature cycles of moderate severity* 
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Piclcett and Grigory [2] reviewed the fatigue evaluation 
procedures and the fracture mechanics approach to fatigue life 

analysis, They made important remarks on crack initiation, 
crack propagation end about the fatigue test specimens. 

They mentioned that specimens loaded in bending are better 
compared with axial or push pull type specimens. They con- 
cluded that the surface strains at the test section on bending 
specimens shake down to stable values after a f ew cycles and 
„ can be measured easily. 

Host studies of low cycle fatigue behaviour have been 
concerned with evaluating the effects of cyclic frequency and 
holding time at peak strain or stress, 

Conway, Berling and Stens £33 did an experimental work 
using servo controlled, hydraulically actuated fatigue machines 
which subjected inductively heated specimens to a push pull 
type of loading. They developed a diametral strain extenso- 
meter and the signal from this device and from the load cell 
in series with the specimen was' fed to a specially designed 
computer network to provide for axial strain control. 

They reported that hold period at peak stress or peak 
strain leads to reduction in the cyclic fatigue life. They 
also reported that the time to fracture is proportional to the 
length of the hold period. This effect is just opposite to the 
number of cycles to . ■ fracture. As the length of the hold 

period increases, number of cycles to cause failure decreases, 
Saraph, Kushwaha and ICakodkar [4] worked on transient 
thermal stress analysis of pressure vessels using finite element- 
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method o They studied the problem of thermal shock, one of 
the most severe loadings experienced by light water cooled 
reactor vessel, when it is subjected to emergency cooling. 

They modelled the vessel as an axisymmetric geometry with 
eight noded isoparametric finite elements. The transient 
temperatures were calculated by Crank-lTicholson' s scheme. 
Detailed temperature and stress plots were given. The 
authors had taken a cylindrical vessel with a hemispherical 
head which is the most or' ■ one in nuclear industry, 

1.2 PRES EM 1 WORK 

In design of vessels subjected to mechanical (pressure) 
and thermal loads, both should be given proper importance. 
Because, the principles upon which thermal stresses can be 
mitigated are many times different to those involved in 
designing for mechanical load, for example, making of 
sections of large thickness can reduce the operating stresses 
for mechanical loads. But this approach may Increase thermal 
stresses, as high thermal gradients are set up in relatively 
thick vessel walls. 

In chemical and petrochemical industries, pressure 
vessels with hemispherical, elliptical, torispherical and 
tor iconic al heads are quite common,' Hence in the present 
work, thermal fatigue analysis of pressure vessels with four 
types of heads mentioned above has been studied using the PEM. 
Thermal stress calculations need prior knowledge of tempera- 
ture distribution. Both temperatures and thermal stresses 
have been determined by the finite element method. Although 
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still it is more popular to use finite difference method for 
temperature calculations but here F3M has been usedo As FEM 
is very straight forward method and implementation of boun- 
dary conditions is an extremely simple affair In addition 
use of FEM for temperature and thermal stress calculations 
reduces the overall data, preparation as same mesh is used 
for botho 

It is known that higher order elements give the same 
accuracy with fewer number of elements. This also requires 
much less data preparation. In addition pressure vessels 
have many crirved surfaces. So the isoparametric eight noded • 
quadrilateral ring elements have been used for temperature 
and thermal stress calculations. 

In the present work, pressure vessels are supported 
by skirts. Insides of the vessels are subjected to cyclic 
temperature variation of heating, ' holding and cooling and 
outsides are insulated. 

Allowable number of cycles as function of heating/ 
cooling rate and hold time are calculated. It is found 
that maximum stress variation occurs in the skirt near the 
head and is the controlling factor for the life of the vessels. 
So design modifications are considered to reduce these stress 
variations in the skirt and thus improve the life of the 
vessel. 
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Finite element method is need .for the thermal fatigue 
analysis of pressure vessels ir. this work. Calculation of 
the thermal stresses requires the temperature dis trihut ion. 
Fii*st the method to determine the temperature distribution 
is given and then the method to determine thermal stresses 
is given. Details of these methods are available in the 
text [<5j but are included hero to make the thesis self- 
contained. 

As the ain of iho thesis j a to study the thermal fatigue 
life of pressure vessels with various heads, design of vessels 
with elliptical, torispherical and toriconical heads is also 
given in this- chapter* 


2.1 QRAITSIEKT TZ -I-khifUIlf) CnhCUh jflPIGll 

Pressure vessels under consideration are axisymmetric 
bodies subjected to axisymmetrie temperature distribution on 
the inside and insulated on the outside* Thus the temperature 
distribution is independent of Q coordinate. Therefore, the 
governing differential equation for the temperature T is 0Q 

i 2T' . i • 2'a 


£ 2 T 


Sr 


■* r ir 4 


a «>t 


O) 


where a is the' thermal dif "uslvlty » 

For this problem f Ini Lc a3 smarts become ring elements. 
To take care of curved botaxd-.rJ.ee, isoparametric eight noded 
quadrilateral , elements arc taken us shown in Fig. 1. 
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temperature distribution over the element is taken as 

= LLJ (2) 

where are the interpolating functions and (T}^ nc ^ are 
the temperatures at the nodes. 

Substitution of Eqn. (2) in Eg.n. (1) gives the residue 
of the element, this residue is minimised by the G-alerkin 
method and finally one gets 

[0] Ce) {ij (no) + [K] (e) {Tj Cne) = {4} (lle) (3) 

■ where 

[0] (e) = ClIiUTJ dr da 

Kl^ = //r({i;,rj L^.rJ + {H,zj IJI.s]) dr dz 

* 

and 

{gJ Uo) = /.{::] (■=£) ran ' (4) 

[C] (e) is the thei’nal capacitance matrix, is the 

thermal stiffness matrix and {qj^ ne ^ is the boundary heat 
flux matrix.- 

2.2 THERMAL STRESSES f m-l K >P J HJhnf ICH 
0?he potential to be minimised is 

l 

I = l-fre dv - /I'Y^dV - Js^ds (5) 

Here the first term represents the strain energy* the second 
term represents the work done by bodj forces and the Ihird 
term represents the’ work done surface forces. 

The pressure vessels under consideration are axisymmetric 
bodies subjected to axisymmetric pressure and temperature loadings. 
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Thus the displacements, strains and stresses are independent 
of 9 coordinate', Thus for this problem also finite elements 
become ring elements as for temperature distribution problem. 
Same isoparametric eight noded quadrilateral elements shown 
in Fig, 1 are taken for this problem, but here each node will 
have the degrees of freedom u and v, 

The displacements over the finite element are taken as 


0 


© © O © 


0 Ij 0 © • © OO0 

1 


^ = D'l (r » a )] ; 


Here u is tbo displacement in the r-direction and v is the 

i 

■displacement in the s-direction, i 

The strain matrix for this axicymmetric problem is » [63 


W- 


Gg 

e~ - _ 


in + 

dz 


Using Eqn, (6), the strain matrix {€} , kqn, (7) becomes 


= [L] {a} (ue) = [B] 1 
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where 


5L i 

0 

1 f.; 

r u 

; i 0 

0 

£ s 

$ 7 ,. 

l 

ah ± 




( 9 ) 


The stress-strain matrices are related to each other as, {jO, 


where 



W = dj]({c{- 

is the elasticity mo tri;: and do the 


( 10 ) 

thermal ’• 


strain matrix, and these are , 

r i 



e(i - « ) 

Tit 5')iterj 


O $ 
T* TT5 


1 


T^3 


0 

o 


sym 1 0 


1-20 

2TP5J 


oc t 




c< >[> 

* T 


cm 


( 12 ) 


l °J 

Here is the coefficient of linear thermal expansion and T 
is the temperature rise * 
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Using Eq.ns# (6), (8) and (10), Eqn. (5) becomes, 

I = i fcil (ne) (2*) fflBf [D] [B] rdrdz {u} (ne) 

- {uj^ ne ) 27T JJl'B'} ~ [D] {€ q } rdrdz 

- |uj (ne) 27* jj^] 1 prds 

- (iij^ ne ^ ZrtjtflX]’* Prdrdz (13) 

By applying the Rayleigh Rits method, one gets, 

[K] (e) fu} (ne) ={l 1 0 }^ ne ' ) + {*' G } (ne) + {^b]^ 316 ^ (H) 


where 

[K] (e ? = znff [B] 1 [D] [3] rdrdz 

{E 0 } (ne) = 27 T^[b] T CDj{€ 0 } rdrdz 

{pj (ne ) = T prdB 

{E b 1 (ne) = 27^11] T Prdrdz (15) 

Here is the stiffness matrix, {Pq}^ 118 ^ is the thermal 

load matrix, consistent surface load matrix and 

{E B } (lle) is body force matrix. 

In this present problem surface loads are the internal 
pressures acting normal to the surface as shown in Pig, 19# 

The consistent surface load matrix is , 

{E} (ne ) =2nJ[N) T |pfj^3 (16) 

If the pressure is acting along the edge ^ = + 1 , then ’ 

_ dz dr 

p r '* p dB d p z “ " p dB 

or p^ dB - pdz 53 P ^ and p^dB = - pdr = - p 
Substituting Eqn# (17)-inEgn# (16), one gets. 



( 17 ) 
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{p}< ne ) 



. 

dr 


r a 





(18) 


Similar expressions for other edges' is ^ = -1 and £ = + 1 
can he written,, Yaltie of this integral is determined by Gauss 
Legendre numerical integration,, 

Body force in the present problem is gravitational 
force. Thus the body force matrix becomes 9 



(19) 


I. 


2.3 


-*“* /*< ^ -n-r> 


2,3,1 Determination of Matrices 

All the matrices developed in Sections 2.1 -and 2.2 are 
evaluated by Gauss Legendre numerical integration scheme. In 
this scheme the quadrilateral is first transformed into a 
square of side 2. Then the function inside the integral is 
evaluated at each Gauss point, and multiplied by the corres- 
ponding weight. Sum of these products give the value of the 
integrals and thus the matrices. 2x2 point integration has 
been used in this work. 

Determination of thermal load matrix {I'q} of Eqn. 

(15) needs some comments. In this matrix, one needs the 
temperature rise T. Reference [6] suggests that this be 
taken as ^average temperature rise of the nodes of the element. 
This does not seem to be an accurate method^ for evaluation of 
this matrix. It is more accurate to determine the temperature 
rise at every Gauss point?; and use it while doing the numerical 
integration. 
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2.3.2 Temperature Determination 

Bq.Ho (3) are first order simultaneous differential 

equations. They are solved here by finite difference 
• iVH 

method. Thus T^ at any time at i node can he 'expressed 

as, 

mP+ 1 _ mP 

m _ _1 1 

x “ At 

and 

T i = 9 T| +1 + (1-9)T? (20) 

1 t 

where p refers to previous time step and p + 1 refers to 
the current time. 

Using Eqn. (20), Eqns (3) become. 


(~x [C] (e) + 9 [K]^ e; ) {T} 


(e) 


,(ne) 


P+1 = [C] (e) - (1^).[K] (e )){T} 

- M 


or 


[e](®) ^T}( ne ^ | P+1 = [L] (e ^{T} (ne ^ 


(ne) jp . 

(21) 

i(ne) ip 

(22) 


These matrices M and [iij can be evaluated once 9 and At 
are fixed. These equations are assembled for the whole domain 
and boundary conditions are applied. The assembled equations 
are simultaneous algebraic equation and are solved by Gauss 
elimination method. This way temperatures at all the time 
steps of the cycle are determined. 

Value, of 9 was talcen as 1.0 in this work and is discussed 


(ne) 


in Chapter III. 
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2.3*3 Thermal Stress Determination 

Finite element eqns* (14) are assembled in the usual 
way and boundary conditions applied* (These are simulta- 
neous algebraic equations and are solved by Gaussian eli-. 
mination for displacements. Stresses are determined using 
Eqn. (10). 


2.4 DESIGN OF VARIOUS PUBS SURF VESSEL HEADS 

A design of the pressure vessel with hemispherical 

heads done by Bharat Heavy Plate & Vessels has been available 

with Department of Mechanical Engineering, Indian Institute 

of Technology, Kanpur. Elliptical, toricpherical and tori- 

conical heads are designed for the vessel using ASME code 

given b;y Chusc & Eber [ll^- 

Design details available for vessel with 
hemispherical heads 


Design pressure 

Working pressure 

Code followed 

Material of the vessel 

Design stress intensity of 
the material 


28 MPa 
26 MPa 

ASME 1980-81, SEC VIII 
SA 336 Gr F11 
1 55 MPa 


9 o 4936x1 O** 6 m 2 -/sec. 


12.474x10' 
7860 kg/m' 


-6 


Yield strength 275 MPa 

Ultimate tensile strength 480 MPa 

Thermal diffusivity 
Coefficient of linear expansion 
Density of the material 
The Fig. 2 gives the other design details of this vessel 
with hemispherical head,. The thickness of the hemispherical 



15 


head has been calculated by the following relation [ll] 


t = 


P R 

2SE - 0.2P 


(25) 


where P = internal pressure 

R = inside radius of head as shown in Pig. 2 ' 

S = allowable stress 
E = joint efficiency. 

Elliptical head: Thickness of the elliptical head is 
• given by [11] 


t = 


PD 

2SE - 0 o 2P 


(24) 


where D is the head skirt diameter. This head has been 
designed for a major axis to minor axis ratio of 2:1. Thick- 
ness comes out to be 140mm. Other details are shown in Pig. 3* 
Torisphorical head : The thickness of the torispherieal • 
head is given by [11] 


_ 0.885 PL 

“ SE - 0.1P 


(25) 


where 1 is the inside crown radius. I lias been taken 974mm 
in this design. With this, head thickness comes out to 155mm. 
Knuckle radius has been taken as 6$ of the inside crown radius 
as required by ASME code [ll] . Other details are given in 
Pig. 4. 


Tor iconic al head : The thickness of the torlconical head 
is given by [ill 

PD- 

( 26 ) 


2cos (SE-0.6P) 

where D-j is the inside ’cone diameter’ at the point 'of f angency 
to the knuckle, ^ is one half of the apex angle. In this 
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design has “been taken as 1435mm and < as 50°. Inside 
knuckle radius has been taken as 6 % of the outside diameter 
of head skirto Other details are given in Pig. 5o 



CHAPTER. Ill 


RESULTS and LI5CUSSIC1T 

As already mentioned, finite element method has been 
used for the study of thermal fatigue life of pressure 
vessels with different types of heads. The stress states 
in various heads have been studied. The effect of rate 
of change of temperature and hold period on the thermal 
fatigue life of the pressure vessels has also been studied. 

5.1 CHECKING THE VALIDITY OF THE PROC-RjiK DEVELOPED 

To check the validity of the finite element program 
developed, a test problem given in £53 was used. Also the 
problem of thick sphere under internal pressure and steady 
state temperature given in GO was used. A good agreement 
of the results • with the available solutions was observed, 

5.2 TRAITS IMT TEMPHRATUR E PIS TR IBUT ION 

Eor transient temperature calculations, both explicit 
(9 = 0) and implicit (9 / 0) were tried with, to find the 
value of 9 which will give minimum or no oscillations of 
temperatures. Generally it is accepted that 9 greater t han 
1/2 Mil give less oscillations £6j. Values of 9 as 0, 

^ (Galerkin’s achemc) and 1 (fully implicit) were tried. 

It was found that fully implicit scheme of 9 = 1 gave 
almost zero oscillations. So in this work this scheme has 
been used through out. 


17 
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3.3 S [TRESS C O MPUTATION iiKD ASSESSFgNT OP FATIGUE LIFE 

After calculating the temperatures in the vessel at 
various time steps of the thermal cycle, stresses were 
determined. [The temperature and stresses in the vessel 
were determined using 2x2 point integration. It was 
observed that the stresses in the vessel reached an extreme 
at the end of the heating portion of the cycle and reached 
another extreme at the end of the cooling portion of the 
cycle. This was found to be true for all types of thermal 
cycles i.e,, with or without any hold period. As only 
maximum alterating stress amplitude was needed for the 
fatigue life evaluation, stx'esses were calculated only at 
these two points of the thermal cycle i.e., at the end of 
heating and cooling portions. In this way, a lot of computer 
time was saved by avoiding the stress calculations at every 
instant of the thermal cycle. 

Further it was observed that the principal stress . 
directions at the above extreme points of thermal cycle were 
different. So to find the stress amplitude, the method 
given by ASME code 02] was used, which is as follow: 

The stresses c - r , c~ z an <3 Xvz a1t: ’ fclle ' fcw0 ex ' treme 

points are calculated and subtracted. These stress diffe- 
" * 
rences are denoted by , c-e'* <rz’ "Exz* Ilie 

cipal stresses corresponding to these stress differences 

are denoted by rj> <r£, <r\ and t:heir differences as 

<r i x - <r z ' > tf- 2 * “ **3* r<md <r 3* “ *V C Half o£ 

maximum difference gives the maximum stress amplitude. The 
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allowable number of fatigue cycles are determined using the 
design 'fatigue curve given in ASMS code £ 12 ]} * 

3*4 STRESS ST-dTE IN VARIOUS VESSELS 

All vessels are subjected to an internal pressure of 
26 MPa* Insides of the vessels are subjected to tempera- 
ture cycles, typical one is shown in the Pig. 6, and out- 
sides are insulated* Temperatures vary from 5°C to 245°C 
on the inside* Two rates of heating/cooling of 1°C/m. and 
1*5°C/m are studied, fur each heating/cooling rate five 
hold periods of 0,1, 2, 3 and 4 hours are considered* The 
thermal cycle with heating/ cooling rate of 1.0°C/min and 2 
hours of hold period is called here as a typical cycle. 

3»4. 1 Stress Distribution hue to Mechanical Loads 

Maximum' principal stress in upper head, cylindrical 
portion and in lower head is shown in Table 2. 

Table 2: Maximum Principal Stress, MPa 


Vessel with 

Upper Head 

Lower Head 

Cylindrical 

Hemispherical 

Head 

152 

173 

136 

Elliptical 

Head 

184 

182 

142 

Torispherical 

Head 

175 ! 

189 

141 

Torieonical 

Head 

117 

126' 

141 
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It is seen that the principal stresses are more in the 

lower head than in the upper head except in the vessel with 

elliptical heads* These stresses are essentially same in 

the cylindrical portion. Principal stress in upper head 

is maximum in the elliptical head and minimum in toriconi- 

cal heado Principal stress in the '.lower head is ma ximum 

in the torispherical head and minimum in toriconical head. 

It is also observed that this maximum principal stress 

occured on the inner side for all the heads. 

Maximum principal stress across AA in the upper heads 

is shown in Pig. To It is seen -that it essentially remains 

constant across the thickness in hemispherical head, and 

fall very sharply in the other heads. 

3.4*2 Temperature Distribution 

Tempera turo distribution in the vessel with spherical 

heads ‘at the end of heati.ij is shown in Pig. 8 and at the end 

of cooling in Pig. 9. Those are for a typical thermal cycle 

in the similar ’./ay temperatures are calculated for other 

vessels and subsequently used for stress calculations. 

3.4.3 Stress Distribution Due to Mechanical and Thermal Loading 

Maximum principal stresses in upper head, cylindrical 

portion, lower head and skirt at the end of heating and 

cooling are given in Table 3« -These are for typical thermal 

✓ 

cycle • 

It is observed that generally maximum principal stress 
occured in lower heads both at the end of heating and cooling. 
Also it is observed that n.- rimum principal stress occured at 
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the end of cooling in all sections of the vessles except in 
the skirt. It is seen that the maximum principal stress in 
the cylindrical portions of all the vessels at the end of 
heating remained more i r k.-' c.anc . TK, same is- observed 
at xh. end of ec\lin 0 , I a m d he .noted that maximum 
values at ifche end of heating and cooling many a times do not 
occur at the same point. 

Maximum principal stress across AA in the upper heads 
is shown in Pigs. 10 and 11. It is observed that maximum 
principal stress in all the heads occuried at the inside 
surface, and falls sharply across the thickness,, 

Maximum principal stress across the cylindrical portion 
BB of the vessel with spherical head is shown in Pig. 12. 

In Figs. 13 and 14 the variation of maximum principal 
stress in vorior.s heads i: a typical thermal cycle is shown. 
(These can be referred to as typical stress cycles produced 
in various heads and are drawn at a point’ in an cl ©lent 
whore stress amj'l.i tv<' c Js observed to be maximum. It is 
observed from these Cijn.*;r that the hemispherical head is 
having the least stress amplitude. It is observed that the 
stress amplitude is increasing in the order of elliptical, 
torispherical and toriconical heads. In Pig. 15 the varia- ■ 
tion of maximum principal, stress in a typical thermal cycle 
in the cylindrical portion of the vessel with hemispherical 
heads is shown. 


# * 
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qtgHL Mm* 
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It is observed from th.se graphs that during the heating 
portion of the cycle, the thermal stresses in the element 
considered are compressive and - during the cooling portion 
they change to tensile stresses . 

3o5 THERMAL FATIGUE LIFE 

In this study, two heating/cooling rates of 1.0 and 
1.5°0/m with five different hold periods as detailed in 
Table 1 are considered. Thermal stresses and maximum stress 
amplitudes are calculated as explained in Section 3.3* 

Maximum stress amplitudes for these thermal cycles 
for vessels with different heads are given in Tables 4 
and 5 for 1 and 1.5°C/m respectively. It is seen that these 
amplitudes hardly change in the upper head and cylindrical 
portions, slightly change in the lower heads and considerably 
change in the skirt regions, for all the vessels. Also it 
is seen that these amplitudes are more in the lower heads 
than those in the upper heads due to the proximity of skirt 
near the lower heads. It is observed that maximum stress 
am plitude always occurs in the skirt region near the vessel* 
Values of these amplitudes are shown in Pig. 16 for •various 
vessels for different hold periods. As expected the maximum 
stress amplitudes increase with increase in hold period and 
heating/cooling rate. It is seen that these amplitudes in- 
crease in the order of hemispherical, elliptical, tori- 
spheric al and toriconical head vessels* 
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Fatigue life of various vessels i.e., allowable number 
of fatigue cycles is calculated using the design fatigue curve 
given'in ASME code jj23« Fatigue lives of various vessels for 
different beating/cooling rates and hold periods are shown in 
Table 6. For heating/cooling rate of 1°C/m these values are 
also given in Fig 0 17 as functions of total cycle time 0 For 
heating/cooling rate of 1 and 1.5°C/m these values are given 
in Fig* 18 as functions of hold period,, For vessels with 
hemispherical heads, with heating/cooling rate of 1.0°C/min, 
the fatigue life decreased by 23$ from a no hold period to a 
hold period of four hours. For the same conditions mentioned 
above, the vessels with elliptical, torispherical and 1 tori- 
conical heads showed a decrease of 28$, 34$ and 36$ respectively. 

With a heating/cooling rate of 1.5°C/min, a change in 
hold period from zero to four hours, vessels with hemi- 
spherical, elliptical, torispherical and torieonieal heads 
show a decrease in fatigue 1 ife of 23$, 29$, 34$ and 46$ 
respectively. 

For a hold period of 2 hours, change of heating/cooling 

« 

rate from 1.0 to 1.5°Q/m vessels with hemispherical, elliptical, 
torispherical and torieonieal heads show a decrease in the 
fatigue life by 50$, 50$, 26$, 14$ respectively. 

3.6 METHODS TO REDUCE THE STRESS AMPLITUDE Iff THE SKIRT REGION 

Throughout the analysis, it is observed that high stress 
amplitudes are developed in the skirt region adjacent to the 
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head. The existing deaden of the skirt caused large thermal 
gradients at tho skirt vessel function, thereby producing 
high stress amplitudes. Thus in the entire analysis the 
skirt is found to be the predominant factor in reducing the 
fatigue life of the vessel. Because of the skirt, higher 
stress amplitudes arc produced in the bottom head also. So 
efforts are made to reduce tho stress amplitude in the skirt, 
thereby increasing the fatigue life of the vessel. The two 
methods which successfully reduced the stress amplitude in 
the skirt are as below: 

1. Changing the position of the skirt. 

2. Providing a suitable temperature distribution 
in the skirt. . 

In the existing designs the skirt is attached at the 
head-cylinder junctions. In the new design, the position of 
the skirt is changed and ' attached only to the head region 
as shown in the Pig. 2. This change in design is studied 
in the vessel with hemispherical head only. For the no hold 
time thermal cycle with 1»0°C/min heating/cooling rate, the 
new design reduces the stress amplitude by about 10$. But 
even this 10$ decrease in stress amplitude does not increase 
the fatigue life much. Co a second method mentioned above 
has been tried with. 

In this second approach, a favourable temperature dis- 
tribution is provided at the skirt to decrease the existing 
high thermal gradients'. Hence one fourth length of skirt 
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adjacent to the head has been subjected to the same thermal 
cycle as that of the inside of the pressure vessel. The 
required temperature variation wi'uh the desired rate of 
change of temperature may be brought by using electric 
heating coils around that skirt. This change considerably 
decreases the thermal gradients at 'idle skirt vessel junction 
and hence the stress amplitude. For a heating/cooling rate 
of 1°C/m, maximum stress amplitudes in upper head, cylinder, 
lower head and skirt of various vessels are shown in Table 7. 
Results for these cases without skirt heating are also given 
in this Table. Maximum stress amplitudes decrease considerably 
in lower heads and skirts. Decrease in skirts is almost 
dr ama tic. Dor vessels with hemispherical, elliptical and 
toriconical heads skirt is no more the -design criterion. 
Fatigue life for vessels with hemispherical and elliptical 
heads become infinite. Fatigue life for vessels with tori- 
spherical and toriconical head increase by 5 0 $ and 45 i» res- 
pectively. It is felt that fatigue life of torispherical 
head can be further improved by shifting the skirt and modi- 
fying the heating rate of the skirt. 



CONCLUSIONS 


In this thesis thcrnal fatigue life of vessels with 
four types of heads namely hemispherical, elliptical, 
torispherical an’d tericonical heads have been studied® 

Insides of the vessels are subjected to heating, holding 

i 

and cooling and outsides insulated® 

It is observed that maximum stress variation occurs 
in the skirt region near the head and is the dominant 
factor for the thermal fatigue life of the vessels with 
all heads. Of the four types of vessels studied, the 
vessel with hemispherical heads has the highest thermal 
fatigue life. The thermal fatigue life of the vessels, 
decreases in the or do: of elliptical, torispherical and 
toriconical heads. 

If stress variation only in the heads is considered, 
it is again observed that hemispherical heads have the 
minimum variation. 

It is seen that thermal fatigue life decreases with 
increase in heating/cooling rate and hold time as expected. 

The decrease is least in the vessel with hemispherical heads. 
The thermal fatigue life decreases in the order of elliptical, 
torisphorioal and toriconical heads. 

Two methods were studied to reduce the variation of 
thermal stresses in the skirt : near the head. One was 
changing the position of the skirt and other was to subject 
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the one fourth of the shirt adjacent to the head to the 
same temperature variation ac that of the inside of the 
vessel. She first method does not reduce stress variation 
much. The second method reduces variation considerably. 
This reduction is dramatic in case of spherical and 
elliptical heads o 
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'able 1 : Thermal Sycles Studied 
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able 3: Maximum Principal Stress, MPa 





























Maximum Stress Amplitudes (MPa) in Various 
Pressure Vessels with heating/cooling rate 






































Sable 5: Maximum Stress Amplitudes (iZPa ) in yarious 
Vessel Heads with a heating and cooling 
rate of 1o5°C/min, 
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Table 6: Maximum Stress Amplitude and Fatigue Life 
c£ pressure Vessels with. Varims Heads 
for Different hold periods. 
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Table 7: Maximum Stress Amplitude (MPa) without and with 
Skirt heating (Heating and Cooling for 4 hours 
at 1 °c/min. , hold period nil, temperature 
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Fig. 3 Pressure vessel with ellipsoidal head. 











Fig. 5 Pressure vessel with toriconical head. 
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Thickness, mm 

Fig. 7 Maximum principal stress across AA in upper heads due to 
mechanical loading. 




Fig. 8 Temperature distribution at the end of heating for 
vessel with hemispherical head. 
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Fig. 9 Temperature distribution at the end of cooling for the 
vessel with hemispherical head. 
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Fig. 10 Maximum principal stress across AA in the 
upper heads due to mechanical and thermo! 
loading. (Typical Thermal Cycle) 
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Fig. 11 Maximum Principal stress across AA in the 
upper heads due to mechanical and thermal 
loading (Typical Thermal Cycle) 
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Fig. 12 Maximum principal stress across BB in the 
cylindrical portion due to mechanical and 
thermal loading (Typical Thermal Cycle) 
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Fig. 13 Stress cycle at a point in the heads giving maximum stress 
amplitude. 
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Fig. 16 Maximum stress amplitude vs hold period 
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